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ENERGY SPECTRA IN SHALLOW U.S. COASTAL WATERS 


by 
Edward F. Thompson 


I. INTRODUCTION 


Wind-generated gravity waves, with periods ranging from 1 to 25 sec- 
onds, are a primary concern in most coastal engineering work. Wind waves 
relate directly to such processes as longshore sand movement, onshore- 
offshore sand movement, and forces exerted on coastal structures. Hence, 
wind waves must be considered in virtually all coastal engineering en- 
deavors, especially in design and maintenance of harbors and exposed 
Navigation channels and in shore protection efforts. Edmisten (1978) 
presents an extensive list of major applications for coastal wind wave 
data. 


Coastal engineers and researchers have wrestled for decades with the 
difficult problem of how to reduce the complexity of ocean wind waves to 
simple, useful approximations. The most common approach has been to rep- 
resent the dominant waves in a particular sea state with three simple 
parameters: Significant wave height, significant wave period, and wave 
direction. The latter parameter is normally unavailable, especially with 
gage data; however, it can be estimated from wind and weather records for 
days of special interest. 


Practical techniques are available for estimating significant height 
and period from a wave record. All of the accepted techniques produce 
reasonable estimates for significant height. However, estimates for 
Significant period are generally more variable and their significance 
is often questionable. One of the best approaches to estimating sig- 
nificant or dominant period is to compute the distribution of wave energy 
as a function of frequency (frequency is the reciprocal of period) and 
use the wave period corresponding to the frequency of greatest energy 
density. This period has a clear dynamic significance. The distribution 
of wave energy as a function of frequency is commonly referred to as the 
wave energy spectrum. 


Although significant height, period, and direction provide a useful 
approximation to many sea states, these simple parameters omit details 
of the sea state which can be important in some applications. For ex- 
ample, the simultaneous occurrence of several prominent wave trains with 
different wave periods and directions is an important consideration in 
relating coastal wave conditions to longshore sediment movement. Design 
of ships and flexible coastal structures, such as floating breakwaters 
and offshore platforms, requires special consideration of how much wave 
energy can be expected at frequencies near the resonant frequencies of 
the structure. Similarly, effective harbor design requires that resonant 
oscillations in the harbor be minimized. Since the effectiveness of 
floating breakwaters in attenuating wave energy depends on the wave fre- 
quency, floating breakwater designs should be formulated and evaluated 
with a knowledge of how wave energy is distributed with frequency. 
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The preceding examples illustrate several applications in which more 
detailed specification of sea state is needed. Much of this additional 
needed detail is provided by the energy spectrum. Evolving technology 
and design techniques in most areas of coastal engineering are generating 
a Similar need for engineers to acquire a more advanced understanding of 
sea State and its climate. 


CERC (and its predecessor the Beach Erosion Board) has operated wave 
gages at numerous U.S. sites during the last 30 years and has acquired 
a volume of shallow-water coastal wave records matched by few other places 
in the world. Most of the gages were mounted on coastal piers in water 
depths of 4 to 6 meters. Significant wave heights and periods from many 
of the gages have been summarized by Thompson (1977). Although none of 
the gages was capable of sensing wave direction, an array of gages was 
used at one site to obtain direction estimates (Esteva, 1977). 


This report presents spectra and spectral summaries for the selected 
CERC gage locations shown in Figure 1; it contains 264 high-energy spec- 
tra which mostly represent shallow-water wave conditions. These spec- 
tra are a significant addition to the existing published literature on 
shallow-water spectra. Field measurements continually show that indi- 
vidual spectra are quite erratic and variable. To identify general 
spectral characteristics, this report contains summaries of approxi- 
mately 7,000 spectra and selected spectral parameters; summaries of 
parameters of the sea-surface elevation distribution function are also 
included. 


The scope of this report is limited. Fewer than one-half of the 
CERC gages for which spectra are available are included. Gages were 
selected for optimum coverage of U.S. coasts and optimum accessibility 
of archived spectra. Only a fraction of the available spectra from 
each gage was considered. For most gages, the spectra represent one 


Presque Isle i 


Figure 1. CERC wave gage locations. 
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reasonably complete year of data. It seems doubtful that more extensive 
summaries would be worth the requisite cost and effort until techniques 
for interpreting and applying spectra are better established. 


This report is intended to provide coastal engineers a better quali- 
tative understanding of shallow-water waves and spectra as well as pro- 
viding quantitative information for engineering use. Qualitatively, it 
provides: 


(a) Perspective on how to interpret shallow-water wave 
energy spectra; 


(b) numerous examples of high-energy, shallow-water field 
wave spectra from the U.S. Atlantic, Pacific, gulf, and Great 
Lakes coasts; 


(c) evidence that the ratio of significant wave height to 
water depth never approaches the commonly used limit of 0.78 
in the water depths considered; and 


(d) perspective on the likelihood of relatively high 
individual waves occurring in succession in a sea state. 


Quantitatively, the report provides the following information helpful to 
coastal engineers: 


(a) Average spectra corresponding to the more common sig- 
nificant height and period combinations at each site; 


(b) frequency of occurrence of sea states with major sec- 
ondary wave trains; 


(c) probability statistics for moments of the sea-surface 
elevation distribution function, which can be related to prob- 
abilities associated with instantaneous surface elevations 
above and below the local mean; and 


(d) a computer program for identifying major peaks and 
valleys in an irregular signal; it has many potential appli- 
cations including identification of major peaks in a spectrum 
and identification of meaningful crests and troughs in a time 
series of sea-surface elevation. 


A brief discussion of the physical characteristics of shallow-water 
ocean waves is presented in Section II. Wave gage data collection, 
analysis, and summarization techniques are described in Section III 
which includes a discussion of the motivation for grouping spectra by 
both significant height and dominant period and computing average spec- 
tra. Spectra computed for several theoretical cases of cnoidal wave 
profiles are presented for perspective in interpreting field spectra. 
The choice, computation, and interpretation of spectral parameters and 
sea-surface elevation distribution function parameters are also dis- 
cussed in Section III. 
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The analyzed and summarized field wave spectra are discussed in 
Section IV which includes samples of spectra representing the highest 
measured wave conditions at each site, average spectra and summaries 
of parameters derived from the spectrum and the sea-surface elevation 
distribution function. Conclusions from the results in Section IV are 
provided in Section V. Major findings are summarized in Section VI. 
Numerous individual high-energy spectral plots are contained in Appen- 
dix A. Average spectra are included in Appendixes B and C. 


II. PHYSICAL CHARACTERISTICS OF SHALLOW-WATER OCEAN WAVES 


Wind waves change considerably when propagating from deep to shallow 
water. Consider a train of uniform waves moving toward shore with crests 
oriented parallel to the bottom contours. The wave height first decreases 
slightly and then increases continually as the waves move into progres- 
Sively shallower water. The phase velocity, group velocity, and wave- 
length all decrease with decreasing water depth. Wave energy is generally 
reduced as waves propagate shoreward. 


It appears contradictory that the wave height in a train of uniform 
waves moving shoreward increases while the energy decreases. Yet the 
two effects occur together mainly because the wave profile changes. The 
profile of a wave in deep water can usually be considered sinusoidal un- 
less the waves are exceptionally high in relation to their wavelength. 
However, the profile of a steep wave in very shallow water is decidedly 
nonsinusoidal (Fig. 2). The troughs are broad and flat, and the crests 
are narrow and high. Deviation of the wave profile from sinusoidal be- 
comes important when the relative water depth (ratio of water depth to 
wavelength) is less than 0.04 for wind waves with very low steepness 
(ratio of wave height to wavelength). Nonsinusoidal profiles also occur 
for relative depths greater than 0.04 when the wave steepness is greater 
than about 0.0015; i.e., for ratios of wave height to the product of 
gravitational acceleration times wave period squared greater than 0.00006 
(U.S. Army, Corps of Engineers, Coastal Engineering Research Center, 1977, 
p. 2-35). For waves with very high steepness (greater than 0.0063 or 
height divided by gT* greater than 0.001), nonsinusoidal profiles occur 
in all water depths. 


Since phase velocity of an ocean wave decreases with decreasing water 
depth, a long wave crest approaching shore at an angle to the bottom con- 
tours has a nonuniform phase velocity along the crest. The velocity dif- 
ferences result in an apparent bending, or refracting, of the crest which 
decreases the angle between the crest and the bottom contours. Refraction 
can affect a wave when the relative water depth is less than one-half. 
Relative depths for the gages in this report range typically from 0.04 to 
1. Hence, refraction has potentially affected many of the wave measure- 
ments considered in this study. Refraction can act to increase or de- 
crease wave heights depending on wave period and direction and bottom 
contour orientation. Since most of the gages are sufficiently nearshore 
that large waves can move bottom sediment seaward of the gage, the bottom 
contour orientation affecting refraction can change with time. 
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Figure 2. Wave profiles and energy spectra for 
several cnoidal wave cases (record 
length: 512 seconds; spectral band- 
width: 0.00977 hertz). 


Another very important phenomenon is wave breaking, which for mono- 
chromatic waves occurs in deep water approximately when the wave steep- 
ness exceeds 0.14 or in shallow water when the ratio of wave height to 
water depth exceeds 0.78. The heights of individual nearshore waves are 
actually quite variable but generally conform to a Rayleigh distribution 
(Thompson, 1974; Goda 1974; U.S. Army, Corps of Engineers, Coastal Engi- 
neering Research Center, 1977). Thus, 13.5 percent of the waves can be 
expected to be higher than the significant height and break farther sea- 
ward than a wave with the significant height. Breaking of a few very 
large waves does not visibly affect the significant height, but when the 
depth-induced breaking starts to affect many large waves, the significant 
height decreases. Quantitative predictions of the decrease in significant 
height as waves begin to break nearshore are given by Seelig and Ahrens 
(in preparation, 1980) for a variety of wave steepnesses and beach slopes. 
Their design curves, based on theory developed by Goda (1975), indicate 
that breaking can begin to affect significant height when the ratio of 
Significant height to water depth is greater than about 0.5. This gen- 
eral observation is consistent with field wave gage data presented by 
Irie (1975). For a typical water depth at the gage of 5 meters, break- 
ing might be expected to have affected all significant heights greater 
than about 2.5 meters. 
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The large effect wave breaking can have on significant height meas- 
ured at most of the CERC gages is shown in Figure 3. The figure shows 
a time history of significant wave height for three Atlantic coast gages 
during a large winter storm. The standard Nags Head, North Carolina, 
gage is represented along with two staff gages on the new CERC Field 
Research Facility (FRF) pier 30 miles north of Nags Head. The FRF gage 
‘farthest seaward shows peak significant heights about 1 meter higher 
than the FRF gage 390 meters farther shoreward. During the peak of the 
storm, Significant heights at the Nags Head gage drop well below heights 
at the seaward FRF gage. The figure suggests that breaking has a major 
influence on significant height at both the nearshore FRF gage and the 
Nags Head gage. 
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Figure 3. Time history of significant wave height for three North 
Carolina gages during a large winter storm. 


For some of the gage sites included in this report, and perhaps for 
all except the accelerometer buoy sites, a sandbar was often seaward of 
the gage. The bar is a somewhat transient feature which can have a major 
influence on the waves by causing high waves to break before they arrive 
at the gages and by acting to reflect energy seaward or trap it between 
the bar and shore. Evidence presented by Wang and Yang (1976) indicates 
the bar acts mainly as an energy dissipator at low tide. At high tide 
for nonbreaking waves, the bar tends to reflect wave energy seaward 
during onshore winds and to trap wave energy between the bar and beach 
during offshore winds. 


Since ocean waves in a sea state are never strictly uniform in either 
height or frequency, a sea state is conveniently characterized by the 
distribution of wave energy as a function of frequency (the spectrum). 
During active wave growth, most energy from the atmosphere is added at 
frequencies slightly higher than the dominant frequency. The energy is 
then redistributed by nonlinear transfers within the wave field to fre- 
quencies slightly lower than the dominant frequency and to very high 
frequencies where it is lost to breaking processes. These phenomena 
result in a gradual shift of the spectral peak to lower frequencies. 
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For waves propagating outside the wave-generation area the nonlinear 
energy transfer to low and very high frequency continues without the 
accompanying replenishment of energy at midfrequencies. Thus, swell 
waves are characterized by narrow spectra with peaks at low frequency. 


There is evidence that the intensity of nonlinear energy transfer 
increases rapidly with decreasing water depth. The nonlinear transfer 
may also narrow the directional spread of energy in shallow-water waves 
(Herterich and Hasselmann, in preparation, 1980). 


Aerial photos of nearshore ocean waters indicate remarkable organi- 
zation in the surface waves (Fig. 4); the photos also indicate that sea 
states dominated by a single visible wave train are relatively unusual. 
Photos typically show two or three apparently independent wave trains 
with different, well-defined directions and wavelengths. The relative 
prominence of the trains often varies with proximity to shore due to 
frequency-dependent shoaling and refraction, and in some cases due to 
redistribution of energy with frequency by nonlinear interactions 
between trains. 


For example, although a long swell wave train often dominates the 
breaker zone in photos taken along the southern California coast, the 
train is invisible two or three wavelengths seaward of the breaker zone. 
The swell waves are preferentially amplified nearshore. McClenan and 
Harris (1975) show several examples of this process occurring. 


The swell waves are also refracted more completely than higher fre- 
quency waves so that at breaking they are often nearly parallel to the 
shore. The water particle velocity associated with low-frequency swell 
waves attenuates slowly with depth. High-frequency waves which have not 
been substantially refracted can travel into the swell-created breaker 
zone with a large component of longshore motion; yet, their particle 
velocities attenuate rapidly with depth. Since longshore sediment move- 
ment requires both a mechanism for initiating sediment movement and a 
longshore drift to transport it, two concurrent wave trains may move 
sediment alongshore more effectively than either train acting alone. 


The presence of a current near a gage site can alter the measured 
characteristics of waves. The longshore current in and near the breaker 
zone affects some of the gage measurements in this report. Tidal cur- 
rents and wind-generated currents may also affect the measurements. The 
Gulf Stream may affect deepwater waves approaching the U.S. Atlantic 
coast by steepening waves enough to cause breaking or, in special cases, 
by refracting the waves so that they never reach the shore (Kenyon, 1971). 
The rate of wave growth in the vicinity of the current may also be 
affected. 


III. WAVE GAGE DATA COLLECTION AND ANALYSIS 


1. Collection. 


The wave gage data used in this report were collected as part of the 
CERC field wave data collection program from 11 gages at’ the 9 locations 
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shown in Figure 1. Table 1 gives some details of the gage locations, 
gage types, and water depth at the gages. The gage designs and charac- 
teristics are described in detail by Thompson (1977), and for the Great 
Lakes accelerometer buoy installations in Thompson (1978). 


The continuous-wire staff gage (manufactured by the Baylor Company, 
Houston, Texas) is the more accurate sensor because the gage directly 
senses surface waves and it interferes very little with waves being 
measured. The continuous-wire gage and the pressure gage (fabricated 
at CERC) give comparable surface wave height estimates when the pressure- 
gage record is compensated for the depth-dependent attenuation of the 
dynamic component of pressure due to surface waves. Suitably compen- 
sated pressure-gage spectra are also comparable to continuous-wire gage 
spectra except at high frequencies (Esteva and Harris, 1970). 


The accelerometer-buoy gage or Waverider (manufactured by Datawell, 
Haarlem, the Netherlands) gives wave heights and spectra which are 
reasonably comparable to those given by the continuous-wire gage for 
frequencies between 0.065 and 0.5 hertz, including virtually all impor- 
tant frequencies in the Great Lakes (Pitt, Driver, and Ewing, 1978). 
The step-resistance staff gage (fabricated at CERC) has been shown to 
have a consistent bias toward high wave heights (Esteva and Harris, 
1970). Significant height estimates from the step-resistance gage may 
be 30 centimeters too high during low wave conditions and 20 percent too 
high during high wave conditions. Aside from the tendency for over- 
estimating wave height and energy, Esteva and Harris show evidence that 
spectra for step-resistance gage records are reasonably consistent with 
continuous-wire gage spectra. 


Time series of sea-surface elevation from staff gages, subsurface 
pressure from pressure gages, and double-integrated sea-surface vertical 
acceleration from buoy gages were transmitted by telephone line to the 
CERC laboratory and used to generate pen-and-ink, strip-chart records 
and digital records on computer-compatible, seven-track magnetic tape. 
Digital records are 20 minutes long. A digital record from a single 
station contains data points at 0.25-second intervals giving a total of 
4,800 data points per location per 20-minute record. Further details of 
CERC field wave data collection procedures are given in Thompson (1974, 
1977) 


Zana lySasts 


a. Spectral Analysis. The digital field wave gage records were 
analyzed by the CERC routine analysis system (Thompson, 1977). One 
record per 6 hours is routinely analyzed; however, one record per 2 
hours was analyzed for some of the Great Lakes data. The routine anal- 
ysis times are chosen to approximately coincide with synoptic times 
(0100, 0700, 1300, and 1900 e.s.t.). The procedure includes editing 
each 20-minute time series to eliminate data points which are obviously 
bad, computing and testing the distribution function of the edited data 
points, and computing the variance or energy spectrum for the time series 
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smoothed with a cosine bell data window. The spectrum is computed with 
a fast Fourier transform (FFT) algorithm applied to 4,096 data points 
(17-minute and 4-second record). 


The spectral analysis procedure assigns a fraction of the total var- 
iance of the record to each of 1,024 frequencies, or frequency lines. 
It is difficult to deal with such finely resolved spectra. To gain sta- 
tistical stability, groups of 11 successive lines in each spectrum are 
combined into bands. The resultant band width is 0.01074 hertz over the 
full range of frequencies considered (0.03 to 1.00 hertz). The uniform 
spectral resolution over the frequency range gives a nonuniform resolu- 
tion in wave period. Significant height is estimated as four times the 
Square root of the total variance assigned to wave frequencies of inter- 
est, usually 0.03 to 1.00 hertz. 


The frequency-dependent responses of some of the gage types required 
some special treatment of the spectra. Accelerometer-buoy gage spectra 
were terminated at the high-frequency end at 0.5 hertz. A special low- 
frequency bound of about 0.065 hertz on the spectrum would have been 
appropriate but was not used because energy at frequencies between that 
and the low-frequency cutoff for other gage types (0.03 hertz) was in- 
Significant in the Great Lakes buoy records. 


Pressure-gage records were compensated for hydrodynamic attenuation 
of the pressure signal with depth by multiplying the energy in each 
spectral band by a frequency-dependent factor. The compensation factor 
becomes very large at high frequencies where the small pressure signal 
is often obscured by noise. Hence, it is necessary to terminate the 
high-frequency end of pressure-gage spectra at a frequency below the 
standard 1l-hertz cutoff for staff gages. The high-frequency cutoff 
used for the Michigan City and Presque Isle pressure gages is 0.33 
hertz. The cutoff for the Pt. Mugu pressure gage is 0.31 hertz. All 
spectral energy assigned to frequencies above the cutoffs was omitted 
from both the spectrum and the significant height estimate. The neg- 
lected energy is unimportant during high wave conditions, but it can 
increase significant height by as much as 0.5 meter during low recorded 
wave conditions. The latter effect is often partially balanced by slight 
overcompensation of the high-frequency energy retained in the spectrum. 


b. Spectral Summarization. 


(1) Averaging. Individual shallow-water ocean wave energy 
spectra are quite irregular and can change substantially and somewhat 
erratically with time (see Fig. 5). Because of these problems, it is 
difficult to identify typical spectra for each wave gage site. However, 
it is desirable to condense the masses of individual spectra into a more 
concise form. 


At this time, no completely satisfactory technique is in use for 
isolating and summarizing general characteristics of field spectra. One 
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Figure 5. 


Energy Density (cm?/Hz) 


00 Of 0.2 0.3 0.4 0.5 
Frequency (Hz) 


Sequence of spectra computed at 6-hour intervals for 
Nags Head, North Carolina, beginning at 0642 e.s.t., 
20 April 1969. H, is 190 centimeters for the top 
spectrum decreasing to 97 centimeters for the bottom 
spectrum. Major spectral peaks are marked with 
asterisks. 
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common approach is to compute spectral statistics which are often based 
on moments of the spectrum, where the nth spectral moment, My LS 
defined as 


By = fs f df (1) 
0 


where f is the frequency and S(f) the spectral energy at frequency f. 
This approach is used in many cases because it is the only tractable way 
to deal with large numbers of spectra. However, some of these spectral 
Statistical parameters have disadvantages which will be discussed later. 
A relatively new and promising approach to spectral parameterization 
involves the use of eigenfunctions (Vincent and Resio, 1977). 


Another approach to the problem of how to characterize and summarize 
ocean wave spectra has been used in conjunction with ship-response pre- 
diction. Hoffman (1974) organized deepwater spectra into 10 groups 
from a site in the North Atlantic according to significant height. For 
each group, the mean and standard deviation of spectral energy in each 
band were computed. The average energy of the one-third highest and 
one-third lowest values in each band was also computed. The importance 
of accounting for variability of spectral shape when predicting ship- 
bending moment responses was demonstrated. 


An averaging procedure for summarizing deepwater spectra from a site 
in the North Atlantic was also used by Gospodnetic and Miles (1974). 
Spectra were grouped according to both significant wave height and 
average wave period. The average spectra were made dimensionless to 
facilitate comparison of spectral shapes. Although most of the varia- 
tions were removed by using the dimensionless format, it was concluded 
that even the dimensionless spectra vary systemmatically with both height 
and period. 


In very shallow water, spectral characteristics would be expected to 
have a pronounced systematic relationship to both significant height and 
dominant spectral period, although the relationship would be somewhat 
obscured when major secondary wave trains exist. Wave height is impor- 
tant because high waves lead to wave breaking in shallow water. Thus, 
any reasonable attempt to compute average spectra at a shallow-water 
location must include stratification of the spectra according to wave 
height or energy. 


The extent of wave shoaling is specified by the relative water depth, 
d/L,, which is directly related to the dominant wave period. Waves in 
very shallow water assume a nonsinusoidal profile which gives rise to 
spectra with peaks at multiples of the dominant frequency. An example 
of the profile and spectrum for several cnoidal wave cases is shown in 
Figure 2; another example is an aerial photo spectrum pair in McClenan 
and Harris (1975, p. 63-64). Because of such systematic effects on the 
shallow-water spectrum, spectra should also be stratified by d/L, or 


by a characteristic wave period. 
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Deviations of the wave profile from sinusoidal occur even in deep 
water for very steep waves as discussed previously. Spectra for such 
waves can also be expected to show energy at multiples of the dominant 
frequency. In such cases the spectral form is related to both wave 
height and period. 


Because of expected systematic influences of wave height and period 
on the shallow-water spectrum, spectra in this study were grouped into 
sets by both significant height and period corresponding to the highest 
spectral peak. Significant height and peak period are convenient because 
they are comparable to height and period parameters used extensively in 
past nonspectral wave analyses. 


Despite the grouping, major differences between individual spectral 
shapes in each set can still occur. Differences between spectra at one 
location for several cases in which the significant height and peak 
period were near the annual means are shown in Figure 6. The differences 
in spectral shape between individual spectra in a set were greater than 
the differences in the mean spectrum between sets in this study and in 
Hoffman's (1974) study. However, any systematic characteristics of each 
set should emerge in the mean spectrum and standard deviation if the 
number of samples is sufficient. This does not imply that the mean 
spectra necessarily represent characteristics. 
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Figure 6. Several individual spectra in the same H, - Tp group 
(Hg = 91 to 122 centimeters; Tp = 8.5 to 9.3 seconds), 
Nags Head, North Carolina. ; 
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Spectra were grouped by 30-centimeter intervals of significant 
height and by variable intervals of peak period (Table 2). Variable 
period intervals were needed because energy is computed for nonuniform 
period bands (but uniform frequency bands) in the CERC spectral analysis 
program. Selection of longer period intervals was necessary so that 
each interval encompasses only one spectral band. If the intervals had 
included more than one spectral band, the averaging procedure would tend 
to produce an unrepresentatively low and broad peak in the average 
spectrun. 


Table 2. Peak period grouping intervals. 
Frequency bands 


Frequency bands Period interval 
| in interval | in interval 
(No. ) 


(No. ) 


Period interval 


1 47 -8 to 8.5 
2 -5 to 9.3 
3 -3 to 10.3 
4 -3 to 11.6 
5 6 to 13.3 
6 3 to 15.5 
6.6 -5 to 18.6 
7.2 6 to 23.3 
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Period intervals of 1 second were used for periods shorter than 6 
seconds (Table 2). Individual spectra with relatively short peak pe- 
riods are generally more poorly focused in frequency than spectra with 
long peak periods. Thus, period intervals encompassing more than one 
spectral band could be used for short periods without creating much 
artificial broadening of the peak in the average spectrum. This treat- 
ment is further defensible on the grounds that (a) short-period waves 
are in relatively deeper water than long-period waves at the measure- 
ment site and hence exhibit less depth-induced modification, and (b) 
spectra with peak periods shorter than 6 seconds are generally not im- 
portant for coastal engineering design applications at ocean sites. 


After the spectra were grouped into height-period intervals, average 
spectra were computed. Intervals were selected for averaging to include 
all intervals containing a reasonable number of cases and all high wave 
intervals regardless of the number of cases. The minimum number of cases 
considered for low and moderate waves was 10 at most locations. The wave 
height above which all intervals were considered varies between locations. 
The average spectrum for a height-period interval was obtained by aver- 
aging energy values in each band. The standard deviation of spectral 
energy values about the mean for each band was also computed. Averages 
and standard deviations of spectra for the locations listed in Table 1 
are discussed in the next section. 


It is important to note that the spectrum for any individual record 
never duplicates the mean spectrum. Any application which is sensitive 
to the precise form of the spectrum should make use of individual spectra 
rather than the mean spectrum. One reasonable approach might be to select 
a small random sample of spectra from each group, apply each spectrum to 
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the problem at hand, and combine the individual results to estimate 
response probabilities, as suggested by Hoffman (1974, 1975). 


(2) Spectral Parameters. 


(a) Number of Peaks. It is relatively common in the ocean 
for two or more independent wave trains with different frequency and 
‘direction to occur simultaneously (e.g., Harris, 1972; McClenan and 
Harris, 1975; Ochi and Hubble, 1976; Thompson, 1977). The marine sur- 
face observation reporting code even includes provisions for reporting 
major secondary wave trains observed visually. Certainly, the concept 
of a simple, single-peaked spectrum is misleading in many ocean wave 
records. Yet, there is a lack of reliable quantitative information as 
to how often major secondary trains occur and how much energy they 
contain. 


To investigate the occurrence of multiple wave trains, a computer 
routine was adapted for use in identifying major spectral peaks. The 
routine (described in App. D) smoothes over spectral peaks and valleys 
which differ in energy density by less than 3 percent of the total 
energy in the spectrum. Examples of individual spectra were shown in 
Figure 5. The major peaks identified in each case by the computer rou- 
tine are marked with asterisks. 


A fraction of the spectral energy is assigned to each major spectral 
peak by a simple method. The spectrum is partitioned at the lowest point 
between successive major peaks (Fig. 7). All spectral energy within the 
partition for a peak is assigned to that peak. The energy is expressed 
as a dimensional quantity and as a fraction of the total spectral energy. 
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Figure 7. Technique for partitioning spectral energy 
and assigning the energy to major peaks. 
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Statistics on the occurrence of multiple spectral peaks identified 
by the above techniques (presented in thé next section) are based only 
on the part of the spectrum between 0.03 and 0.5 hertz for the staff and 
accelerometer-buoy gages. The high-frequency cutoff for the pressure 
gages was 0.33 hertz in the Great Lakes and 0.31 hertz at Pt. Mugu as 
discussed earlier. 


(b) Spectral Peakedness. It is highly desirable to param- 
eterize a spectrum. Various parameters have been considered in the lit- 


erature, but most are derived from moments of the spectrum as defined in 
equation (1). 


In practical application equation (1) becomes 


f 
My, = far S(£) f"d£ (2) 


fon 


where f77 is the low-frequency spectral cutoff, and fyr the high- 
frequency spectral cutoff used for computation. When applied to field 
records of finite length, the integral in equation (2) becomes a 
summation. 


Equation (2) shows that the higher the value of fyp, the greater 
the magnitude of m, provided that S(f) is greater than zero at high 
frequencies. In practice, the choice of fyp has a significant influ- 
ence on the values computed for m, especially for the higher moments. 
The choice of f;, is less ambiguous in practice because contributions 
to m, from the low-energy and low-frequency part of the spectrum are 
generally negligible. 


The sensitivity of m, and parameters based on m, to the choice 
of fyr was investigated by Rye (1977). He considered two theoretical 
spectral shapes: a sharply peaked JONSWAP spectrum and a flat JONSWAP 
spectrum equivalent to a Pierson-Moskowitz spectrum. Rye recommended 
three spectral parameters for general use: 

(a) Significant wave height, H, where 

Hg = 4¥m (3) 

(b) period of the spectral peak, Tp 


(c) spectral-peakedness parameter, ® where 


i) 
Ex 


Q = <> | f[s(#)12aF (4) 
0 


The spectral-peakedness parameter was originally proposed and shown to 
be directly related to the average number of high waves oceurring in 


AT 


succession in simulated wave records by Goda (1970). However, Goda (1976) 
could not find a clear relationship between Q and the extent of group- 
ing of high waves in field records. 


The value of Qp as a spectral parameter is shown in Figure 8. The 
two spectra shown have nearly the same significant height and peak period. 
The spectra also have comparable values of the frequently used spectral- 
width parameter, €, defined as 


2 
m 
eee loa! (5) 
Mo Mm, 


However, Q) for one spectrum is 66 percent greater than for the other, 
indicating substantial differences in spectral shape. 


Energy. 


Frequency 


Figure 8. Comparison of two measured spectra 
from the North Atlantic Ocean with 
Hg = 3.3 meters, Tp = 10.5 seconds, 
and « = 0.61 (after Hoffman, 1974). 
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Both Q and e are designed to fit the concept of a single-peaked 
spectrum, but multipeaked ocean wave spectra are commonly observed. The 
suitability of spectral parameters for representing multipeaked field 
wave spectra has received little attention in the literature, although 
Ochi and Hubble (1976) have developed a computationally complicated tech- 
nique for parameterizing field spectra with two major peaks. 


Since Qp seems to be a potentially important and useful parameter, 
it was computed for each of the spectra considered in this report. 
was based on the part of the spectrum between 0.03 and 0.2 hertz for 
staff and accelerometer-buoy gages, between 0.03 and 0.33 hertz for the 
Great Lakes pressure gages, and between 0.03 and 0.31 hertz for the Pt. 
Mugu pressure gage. 


c. Parameters of Distribution Function of Sea-Surface Elevation. 
The distribution function for instantaneous sea-surface elevations pro- 
vides useful insight on shallow-water spectra and spectral parameters. 
It also provides probabilities associated with instantaneous surface 
elevations above the mean. 


The distribution function of sea-surface elevations can conveniently 
be parameterized by its moments defined as: 


N 
= 3 iM? oh. 6 
qd, sep ie Gg) (6) 
where 
Gn = nth moment of the distribution function of sea-surface 
elevations 
N = number of intervals in the distribution function 
n; = sea-surface elevation associated with the ith interval 
in the distribution function 
p(n,) = probability associated with n,; 
L 


The zeroth and first moments, Gg and q,, are equivalent to the mean 
and variance of the distribution function. q3 and q, are often re- 


ferred to as the skewness and kurtosis of the distribution function. 


The distribution function of sea-surface elevations is often assumed 
to be Gaussian. When normalized, the Gaussian distribution function has 
a mean of zero, a variance of one, a skewness of zero, and a kurtosis of 
three. 


Steep waves in shallow water assume a decidedly nonsinusoidal pro- 
file with broad, flat troughs and narrow, high crests. The distribution 
of sea-surface elevations measured at a point is obviously non-Gaussian. 
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The skewness and perhaps the kurtosis of surface elevations in very 
shallow water can be expected to differ from the values for a normal- 
ized Gaussian distribution. For the cnoidal wave profile in Figure 

9, the skewness and kurtosis are equal to 2.0 and 5.8, respectively. 
The figure also shows the distribution function, skewness, and kurtosis 
for several idealized wave profiles. 
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Figure 9. Distribution function of sea- 
surface elevations for several 
idealized wave profiles. 


Both skewness and kurtosis for the normalized distribution function 
of sea-surface elevations were computed for each record and are summa- 
rized in the following section. For the accelerometer-buoy data, the 
distribution functions represent doubly integrated vertical-surface 
accelerations rather than directly measured sea-surface elevations. 

For the pressure-gage data, the distribution functions represent sub- 
surface pressure fluctuations. 
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IV. SPECTRA, SPECTRAL SUMMARIES, AND SEA-SURFACE 
ELEVATION DISTRIBUTION FUNCTION PARAMETERS 


1. High-Energy Wave Spectra. 


Since high-energy spectra are of greatest concern in most engineering 
applications, plots of individual spectra containing the most energy for 
each station are provided in Appendix A. Spectra are grouped by gage 
location in the order listed in Table 1. The appendix includes 24 spec- 
tra for each location, arranged and numbered in descending order of sig- 
nificant height. 


Each plot in Appendix A shows the significant wave height, period 
corresponding to the highest spectral peak, spectral-peakedness parameter 
(Qo), and number of major spectral peaks identified by the computer rou- 
tine SMOOTH (App. D). Spectral plots for days associated with tropical 
storms and hurricanes are labeled accordingly. The energy density scale 
can vary between plots. 


Pressure-gage spectra in Appendix A are cut off at high frequencies 
above 0.33 hertz for the Great Lakes gages and 0.31 hertz for the Pt. 
Mugu gage. Several of the Pt. Mugu high-energy spectra were obviously 
overcompensated for attenuation of pressure with depth so that the high- 
frequency energy was much greater than the low-frequency energy. These 
spectra were omitted from Appendix A. This problem can arise because 
the high-frequency cutoff is constant while the compensation increases 
as the depth of water over the gage increases with the tide. Thus, the 
high-frequency end of the spectrum is sometimes overcompensated when 
water levels are high. This problem did not arise with pressure-gage 
spectra in the Great Lakes (where there is no tide). 


Energy computed for wave periods longer than 30 seconds (frequencies 
less than 0.043 hertz) does not directly represent wind-generated waves. 
In some cases it is questionable whether such very long-period energy 
has any physical meaning. Therefore, it was omitted from spectral plots 
for the Great Lakes and Pacific coast locations and should be ignored in 
plots for the Atlantic and gulf coast locations. 


Spectral plots for each location in Appendix A are followed by a table 
which contains auxiliary data for each plot, including sequential number, 
date, time, significant wave height, peak wave period, water depth, rela- 
tive water depth (where is the shallow-water wavelength for a wave 
with frequency equal to the frequency of maximum spectral energy density), 
ratio of significant wave height to water depth, wave steepness, ratio of 
depth to gTf (which is related to wavelength), spectral-peakedness pa- 
rameter, and skewness and kurtosis of the distribution function of sea- 
surface elevations. Wavelength used in the tables was computed in each 
case from period corresponding to the highest spectral peak and the equa- 
tions of linear wave theory. The largest ratio of significant wave height 
to depth is 0.54. None of the significant height-to-depth ratios approach 
the commonly used limit of 0.78. 
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Multiple major spectral peaks are evident in many of the spectra in 
Appendix A. However, some of the multipeaked spectra and many of the 
single-peaked spectra appear to approximate the form of the JONSWAP spec- 
trum for waves undergoing active generation (Hasselmann, et al., 1973). 
Only seven of the spectra selected for inclusion in Appendix A were col- 
lected during hurricanes. This small sample does not illustrate any 
clear differences between spectra generated by hurricanes and spectra 
generated by extratropical storms. 


2. Mean Spectra. 
a. Grouped by Height and Period. Spectra from each gage site were 


grouped into sets by significant height and peak spectral period. Inter- 
vals of 30.5 centimeters were used for significant height; intervals of 
about 1 second were used for peak period (Table 2). 


Mean spectra were plotted for all height-period intervals containing 
more than one case with heights greater than a cutoff height, which 
varied with location. In addition, mean spectral plots were generated 
for height-period intervals with heights below the cutoff if there were 
a reasonably large number of cases in the interval. This approach pro- 
vided mean spectral plots for the most statistically meaningful low wave 
cases and for all high wave cases. Energy computed for periods longer 
than 30 seconds was omitted from plots for Great Lakes and Pacific coast 
locations as discussed earlier. 


The mean spectral plots grouped by location are given in Appendix B. 
The approximate relative water depth and ratio of significant wave height 
to water depth for each location are listed in Table B-1; the dates 
covered are in Table 1. 


The range of significant wave height and peak spectral period repre- 
sented and the number of cases are listed on each plot. Each plot repre- 
senting more than one case contains several dashlines in addition to the 
solid line showing the mean spectrum. The dashline nearest the mean and 
above it represents the mean plus one standard deviation. Similarly, 
the dashline nearest the mean and below it represents the mean minus one 
standard deviation. When the mean minus one standard deviation is nega- 
tive, it is plotted as zero. For plots representing more than two cases, 
another dashline above the mean shows the highest single-energy density 
value occurring in each band. This line forms an envelope inside which 
all of the individual spectra would fall. The dashlines should not be 
treated as a spectrum, and are included only as indicators of the observed 
variability of energy density in each spectral band. Some of the plots 
show a relatively large concentration of energy at the lowest frequency 
(0.005 hertz). This concentration is at a frequency which is irrelevant 
to wind-generated gravity waves and should be ignored. 


The mean spectral plots for each location in Appendix B are arranged 


in order of increasing significant height and peak period in most cases; 
i.e., all mean spectra in a colum represent the same peak-period interval 
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and all mean spectra in a row represent the same significant height 
interval. Thus, the effect of increasing height on spectra with similar 
peak periods and the effect of increasing peak period on spectra with 
Similar heights can be easily followed. A small number of plotted mean 
spectra are not included in Appendix B because the spectra (a) usually 
represented low or moderate wave heights with relatively few cases or 

high wave heights with only one case for which the spectrum is shown in 
Appendix A, or (b) were high-energy cases for Pt. Mugu which had obviously 
been overcompensated for attenuation of pressure with depth. 


Mean spectra for low wave heights and peak periods less than 5 seconds 
along the Atlantic coast show evidence of secondary peaks at low frequency 
in many cases. The low-frequency peaks are broad and flat, are generally 
located near 0.1 hertz, and are likely to represent low-energy swell in 
the Atlantic Ocean. Corresponding broad, low-frequency secondary peaks 
are also evident in mean spectra for Pt. Mugu on the Pacific coast. These 
peaks are not evident in mean spectra for the gulf and Great Lakes gages, 
although the Presque Isle pressure-gage spectra show a sharp secondary 
peak which is consistently centered on 0.038 hertz. This low-frequency 
peak at Presque Isle may not be physically meaningful. For all locations, 
the differences between mean spectra are smaller than the differences 
between the individual spectra used to compute the means. 


Mean spectra from the ocean gages for peak periods longer than about 
10 seconds show evidence of secondary peaks at high frequency, especially 
for high wave heights. For example, several selected mean spectra for 
high wave, long-period groups are plotted in Figure 10 in dimensionless 
form. Energy density was made dimensionless for each mean spectrum by 
dividing by the total energy in the spectrum and converting to percent; 
frequency was made dimensionless by dividing by the peak frequency. 


The figure clearly shows that the main secondary peak in each mean 
spectrum is at twice the peak frequency. The spectrum for one cnoidal 
wave (see Fig. 2) is also shown to clearly resemblé the field spectra. 
A section of the Lake Worth pen-and-ink strip-chart record taken at 
nearly the same time as the Lake Worth spectrum in Figure 10 is shown 
in Figure 11. The profiles of the larger waves closely resemble the 
cnoidal wave profile. Many of the high-frequency secondary peaks in 
mean spectra for high, long-period waves are probably indicative of 
cnoidal-type wave profiles. Spectra for high, long-period waves in 
which secondary peaks at about two or three times the peak frequency 
represent independent wave trains are expected to be rare. 


The plots in Appendix B also show interesting variations with wave 
height in mean spectra with short peak periods. In some cases, the 
mean spectra for high waves exhibit secondary peaks at low frequency. 
The emergence of the low-frequency secondary peaks for Nags Head is 
shown by dimensionless spectra plotted in Figure 12. The figure also 
indicates that the largest secondary peak is at one-half the frequency 
of the dominant frequency. A section of pen-and-ink strip-chart record 
corresponding to two high wave, short-period cases is shown in Figure 13. 
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Figure 10. Dimensionless mean spectra for several high wave, long-period 
groups and for the cnodial wave shown in Figure 2, case 3. 
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Figure 11. Section of pen-and-ink strip-chart record from Lake Worth, 
Florida, beginning at 0810 e.s.t., 20 Feoruary 1909. 
Parameters are derived from analysis of a digital record 
beginning at 0723, e.s.t., 20 February 1969. Individual 
waves with height greater than H, are marked by asterisks. 
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Figure 12. Dimensionless mean spectra for several high wave, 
short-period groups at Nags Head, North Carolina. 
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Figure 13. Two pen-and-ink strip-chart records and corresponding spectra from 
Nags Head, North Carolina. (a) Record at 0800 e.s.t., 20 April 1969; 
analysis at 0642 e.s.t., 20 April 1969; (b) record at 2000 e.s.t., 
1 March 1969; analysis at 2040 e.s.t., 1 March 1969.- of 
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In one case the spectrum has a single dominant peak; in the other case 
the spectrum has a large peak at one-half the dominant frequency and 
another major peak at about 1.5 times the dominant frequency. 


A laboratory parallel of this behavior has been observed by Sawaragi 
and Iwata (1976). A single-peaked spectrum was generated in shallow 
water of uniform depth. As the waves propagated and broke, an increas- 
ing fraction of the energy shifted to frequencies which are multiples of 
the dominant frequency, fp. Eventually, the spectral peak at 2fp was 
higher than the peak at fp- 


The individual spectral points used in computing mean spectra are 
plotted for Lake Worth in Appendix B. The distribution of the spectral 
points about the mean for each band is decidedly unsymmetric in most 
cases. The distribution usually has a high positive skewness indicating 
some relatively high values above the mean in the distribution which are 
not offset by correspondingly low values below the mean. However, the 
distribution of spectral points about the mean for the band of maximum 
energy density tends to be more symmetric and more nearly Gaussian. 


Third and fourth moments of the normalized distribution of spectral 
points in selected bands were computed for all locations considered in 
this report. The moments are defined analogous to the moments of the 
distribution function of sea-surface elevations (eq. 6). For all loca- 
tions, both the third and fourth moments are generally lowest in the 
vicinity of the spectral peak indicating a relatively symmetric and 
broad distribution about the mean. 


Typical variations of the third and fourth moments of the normalized 
distribution of spectral points as a function of frequency are shown in 
Figures 14 and 15 for mean spectra with peak periods between 7.2 and 7.8 
seconds at Nags Head. The third and fourth moments in the band of peak 
energy density actually fall below the values for a normalized Gaussian 
distribution. Deviations from the normalized Gaussian distribution over 
all frequency bands are generally less for high-energy spectra than low- 
energy spectra. It may be concluded that the normalized Gaussian dis- 
tribution gives a reasonable representation of the scatter of individual 
spectral energy densities about the peak of each mean spectrum in Appen- 
dix B, but it is not reliable for representing scatter about lower energy 
bands away from the peak. 


Seasonal variations in the mean spectrum for each height-period in- 
terval were studied for Nags Head. Four 3-month seasons were selected 
using the seasonal summaries in Thompson (1977) as a guide. Mean spectra 
were computed for winter (December, January, and February) and for summer 
(June, July, and August). Seasonal differences in the mean spectra did 
not appear to be significant. The number of cases per height-period 
interval was considerably reduced when the spectra were grouped by sea- 
son. Hence, extensive seasonal summarization was not considered in this 
Study. 
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Comparison of mean spectra in Appendix B for the same height-period 
intervals at different locations provides additional perspective. Mean 
spectra for the Atlantic coast locations are generally quite similar to 
each other (Fig. 16,a and b). Mean spectra for the gulf coast location 
differ from the Atlantic coast spectra in that spectral energy at fre- 
quencies lower than 0.1 hertz is absent in virtually all cases. Mean 
Spectra for the Great Lakes locations show an even more pronounced 
absence of low-frequency energy, especially for low wave conditions. 
This effect is illustrated in Figure 16a. The absence of low-frequency 
energy gives the gulf and Great Lakes mean spectra a tendency to have 
higher peaks than the mean spectra for corresponding height-period 
groups at Atlantic coast locations. 


Mean spectra for the Pacific coast locations show that peak periods 
longer than 10 seconds are more common than for Atlantic coast locations. 
Reduced amounts of energy are also shown at frequencies higher than about 
0.15 hertz in comparison to Atlantic coast spectra with long peak periods 
(Fig. 16b). Hence, the Pacific coast long-period mean spectra generally 
have higher peaks than their counterparts for the Atlantic coast. The 
few Pacific coast mean spectra shown with peak periods shorter than 9 
seconds also show a distinct swell energy concentration at frequencies 
less than 0.1 hertz. 


b. Grouped by Height, Period, and Water Level. Numerous field 


studies have indicated that shallow-water spectral shape is dependent on 
local water depth (e.g., Barnett 1969; Hasselmann, et al., 1973; Wang and 
Yang, 1976). Therefore, many of the results in this report are expected 
to be strongly influenced by the water depth at the gage site. The higher 
wave conditions at many of the sites can represent situations where waves 
occasionally break seaward of the gage, or, in the extreme, cases where 
the gage is actually in the surf zone as discussed in Section II. 


To give more perspective on the systematic influence of water depth 
on spectral shape, the spectra within each significant height-peak period 
interval for the gage at Nags Head were further grouped according to mean 
water level in the record. Three water level groups were arbitrarily 
defined. Plots of average spectra by water level for height-period in- 
tervals in which there were a reasonably large number of cases are shown 
in Appendix C. Each plot is labeled with the height and period intervals 
represented and with a mean water level designator. A mean equal to 1 
indicates the mean water level was relatively low; a mean of 2 indicates 
midtide water levels and 3 high water levels. Water level 2 encompassed 
a range of 49 centimeters. The differences between mean spectra as a 
function of mean water level are small. 


Sh Spectral and Sea-Surface Elevation Distribution Function Parameters. 


a. Secondary Spectral Peaks. The computerized procedure described in 
Section III and Appendix D was used to identify major peaks in all spec- 
tra considered for this report. For five of the six ocean locations only 
about one-third of the shallow-water spectra are single peaked (Fig. 17). 
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Figure 17. Percentage of spectra versus 
number of major spectral peaks. 


At the remaining ocean location (Pt. Mugu) only 5 percent of the spectra 
were single peaked. The infrequency of single-peaked spectra at Pt. Mugu 
may be attributed both to common simultaneous occurrence of sea and swell 
waves and to pressure-compensation effects on the spectrum which often 
make the high-frequency end of the spectrum look like a meaningful peak. 
At all ocean locations, about one-half the spectra are double peaked with 
the rest containing three, four, and in a few cases, five major peaks. 
Several spectra with five major peaks and corresponding pen-and-ink 
traces are shown in Figure 18. 


Single-peaked spectra are more common for the gulf and Great Lakes 
sites than for the ocean sites (see Fig. 17). About 40 percent of the 
Naples spectra and about one-half of the Great Lakes spectra from accel- 
erometer buoys 6 kilometers offshore were single peaked. Spectra from 
the Great Lakes nearshore pressure gages show that only 20 to 30 percent 
of the cases are single peaked. However, these percentages are consid- 
ered low because the pressure-compensation procedure often creates a 
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Pen-and-ink strip-chart records and corresponding spectra 
for some cases with five major spectral peaks: (a) Nags 
Head record at 0400 e.s.t., 25 December 1968; analysis 

at 0643 e.s.t., 25 December 1968; (b) Nags Head record 

at 0800 e.s.t., 24 February 1969; analysis at 0642 e.s.t., 
24 February 1969; (c) Huntington Beach record at 1900 
€.s.t., 9 May 1972; analysis at 1840 e.s.t., 9 May 1972. 


small peak at high frequency. Small peaks at frequencies less than 0.05 
hertz appeared in some of the Great Lakes spectra, especially spectra 
from the Presque Isle pressure gages (App. B). These peaks were consid- 
ered spurious and are not counted as major spectral peaks. Double-peaked 
spectra are shown to be common at both the gulf and Great Lakes locations. 


It was noted that secondary spectral peaks do not always represent 
independent wave trains. When waves are nonlinear, i.e., when H,/d or 
Hg/Lp are large, the wave profile will be nonsinusoidal. For such cases, 
a secondary spectral peak at twice the dominant frequency is expected 
even though a photo would show a single train of waves. The effect is 
most evident for high, long waves at low tide; i.e., cases with high 
values of Hsg/d and low values of d/Lp. Appendixes A and B clearly 
show that secondary peaks tend to occur systematically for high, long 
shallow-water waves. However, the joint distribution table of number 
of major spectral peaks versus significant height at each location shows 
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no clear evidence of a relationship (see Table 3). The correlation be- 
tween number of peaks and Hs is less than 0.3 at all locations (Table 
4). The joint distribution table of number of major spectral peaks ver- 
sus period corresponding to the highest spectral peak at each location 
also shows little evidence of a relationship (Table 5). The correlation 
between number of peaks and Tp is also very low (Table 4). For most 
wave conditions the occurrence of secondary spectral peaks does not 
appear to be systematically related to Hg and Tp. The overall sta- 
tistics on multiple peaks as indicators of multiple wave trains should 
be reasonably correct. 


Several other attempts have been made to quantify the likelihood 
of having multiple wave trains. Shepard and Inman (1951) report visual 
wave observations to identify secondary and dominant wave trains from 
the Scripps pier in southern California during a 10-month period. They 
report secondary wave trains in 56 percent of 302 observations. This 
figure may be an underestimate because of the difficulty of visually 
identifying secondary trains, but is reasonably consistent with the 
Huntington Beach gage data. 


Table 3: Percentage of cases in each significant height group 
corresponding to number of major spectral peaks. 


Major spectral pouks 


a) ca a SS eS A] 


Nags Head, N.C. 


INo observed cases for combinations of H, and number of peaks. 
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Table 4. Linear correlation betwecn number of major spectral peaks, 
significant height and period corresponding to the highest peak. 


Correlation 
vs. 


Location Hg Tp vs. 


No. of peaks No. of peaks 


Atlantic City, N.J. 0.12 0.07 
Virginia Beach, Va. 0.10 0.10 
Nags Head, N.C. 0.16 0.06 
Lake Worth, Fla. 0.02 0.05 
Naples, Fla. 0.17 0.05 
Michigan City, Ind. (buoy) H 0.24 0.17 
Michigan City, Ind. (pressure) 0.26 0.24 
Presque Isle, Pa. (buoy) 0.09 0.18 
Presque Isle, Pa. (pressure) 0.19 0.32 
Huntington Beach, Calif. 0.01 0.17 
Pt. Mugu, Calif. 0.02 0.16 


— z —— 


Table 5. Percentage of cases in each peak period group 
corresponding to number of major spectral peaks. 
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Hoffman (1978) reports on the occurrence of multipeaked spectra from 
hindcasts generated for a grid point in the North Atlantic Ocean by the 
Spectral Ocean Wave Model of the U.S. Navy Fleet Numerical Weather Cen- 
tral, Monterey, California. The criterion for accepting a secondary 
peak was that its peak energy be at least 30 percent of the energy in 
the highest peak and that it be reasonably separate in frequency from 
the main peak. About 25 percent of 572 cases considered were multi- 
peaked. There was a tendency for multipeaked spectra to be more common 
during low wave conditions than during high wave conditions. 


It is of some engineering interest to know whether a spectrum with 
a certain Tp value is more likely to contain major secondary peaks 
than a spectrum with a different Tp value. Evidence at many locations 
shows that spectra associated with certain ranges of Tp do indeed ex- 
hibit a preference for multiple peaks. At Atlantic City, Virginia Beach, 
and Nags Head on the Atlantic coast and at Huntington Beach on the Pacific 
coast, multipeaked spectra are more common for long and short T,, cases 
than for intermediate T, cases (shown in Fig. 19). Multipeaked spectra 
are least common for T, values between about 8 and 10 seconds at the 
three Atlantic coast locations. At Huntington Beach multipeaked spectra 
are least common for Ty values of 12 to 13 seconds. At Lake Worth and 
Naples, there is no clear evidence of a relationship between T,, and 
the occurrence of secondary spectral peaks. Data on secondary peaks from 
Pt. Mugu and the Great Lakes locations were not considered because of | 
difficulties with the compensated pressure spectrum and because of the 
small number of Great Lakes cases. 


Additional insight on the occurrence of secondary spectral peaks can 
be obtained from a detailed look at their characteristics for a repre- 
sentative Atlantic coast location and the single Pacific coast location. 
For example, Figure 19 shows that 70 percent of the Nags Head spectra 
with T between 6 and 7 seconds have major secondary peaks; Figure 20(a) 
shows in a histogram that the periods of those secondary peaks range from 
2 to 17 seconds. The figure also shows that secondary peak periods are 
fairly evenly distributed over the range of 3 to 17 seconds. The rela- 
tively low percentage of secondary peaks at 7 to 8 seconds shown in Fig- 
ure 20(a) is artificially low because the CERC analysis procedure gives a 
nonuniform resolution in wave period which tends to decrease the number 
of cases in the 7- to 8-second intervals and increase the number of cases 
in the 6- to 7- and 8- to 9-second intervals. The procedure also precludes 
any ‘cases wath intervals) of 11 tol2, 15) to, 14) 15) to lo. “anid pili7stomzo0 ; 
seconds. 


Although the histograms in Figures 20 and 21 show how often secondary 
peaks occur, no indication is given of how much energy is contained in 
the secondary peaks. Therefore, another curve was added to the figures 
to show the average energy for secondary spectral peaks in each period 
interval. The energy associated with each peak was estimated by the pro- 
cedures discussed in Section III and shown in Figure 7. Relatively high- 
energy secondary peaks tend to occur in Figure 20(a) at periods, slightly 
less than T. and at periods of 12 to 14 seconds. The long-period sec- 
ondary peaks generally represent swell waves. 
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Figure 20. Frequency of occurrence and average energy for 
secondary peaks as a function of secondary peak 
period at Nags Head, North Carolina. 
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Figure 21. Frequency of occurrence and average energy for 
secondary peaks as a function of secondary peak 
period at Huntington Beach, California. 


47 


bh 
oO 
(o} 


200 


Average Energy (cm?) 


Figure 20(b) shows characteristics of secondary spectral peaks at Nags 
Head when period of the dominant peak is 12 to 13 seconds. Secondary 
peak periods of 4 to 9 seconds are common. In most cases these secondary 
peaks represent locally generated wave energy or swell waves generated 
near the Nags Head coast. The average energy associated with these peaks 
is relatively high. Secondary peaks at periods longer than 10 seconds 
are almost totally absent. 


At Huntington Beach, spectra with relatively short Tp have a very 
strong tendency for a secondary peak at long period (see Fig. 21,a). 
Locally generated waves at Huntington Beach are usually superimposed on 
long-period swell waves. Huntington Beach spectra with long Tp have 
a small tendency for secondary peaks with short periods (see Fig. 21,b). 
They also have a small tendency for secondary peaks at long period, which 
presumably correspond to secondary swell waves superimposed on the domi- 
nant swell. However, most Huntington Beach spectra with long Tp appear 
to be dominated by the main spectral peak. 


b. Spectral Peakedness. The spectral-peakedness parameter, > 
defined by equation (4) was computed for each spectrum. The spectral- 
peakedness parameter is shown on each plot in Appendix A. 


Cumulative distribution functions of Q) are shown in Figure 22 
for the selected locations. The Q») values cover a wide range at most 
locations—from less than one to greater than five. values tend to 
be smallest for the Atlantic coast sites. The highest Qp values tend 
to occur for the Great Lakes gages, with the pressure gages having higher 
Q> values than the buoy gages. Q values for the gulf and Pacific coast 
gages tend to be intermediate to values for the Great Lakes and Atlantic 
coast gages except for the highest 5 percent of the Pt. Mugu Qp values. 
The cumulative distribution function for Pt. Mugu shows a rapid increase 
in Qp» values associated with probabilities of less than 5 percent. 
This feature is a result of occasional large overcompensation of the 
pressure spectrum which creates a high spectral peak at high frequency. 
Several cases in which Qp was greater than nine are omitted from Fig- 
ure 22. These cases all correspond to very low wave conditions except 
for several overcompensated Pt. Mugu cases. 


Goda (1976) indicated that high values of are associated 
with strong grouping of high waves. Pen-and-ink record traces for sev- 
eral cases with reasonably high waves and high values are shown in 


Figure 23. Some grouping of the high waves is evident in the records. 
Records with low Qp) values in Figures 13 and 18 show little evidence 
of wave grouping; yet the record in Figure 11 for high, long waves shows 
evidence of high waves occurring together, but the Q> value is only 
Sore 


The spectral-peakedness parameter might be expected to weakly corre- 
late with the number of major spectral peaks so that single-peaked spec- 
tra tend to have high peakedness. The data indicate a low correlation 
at some locations, but not at others (Table 6). Qp) also appears to be 
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Location No. of Cases 
Atlantic City 1037 
Virginia Beach 982 
Nags Head 945 
Lake Worth 982 
Naples 866 
Michigan City (Buoy) 246 
Michigan City (Pressure) 127 
Presque Isle (Buoy) 110 
Presque Isle (Pressure) 119 
Huntington Beach 746 
Pt. Mugu 799 
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Figure 22. Cumulative probability distribution functions 
for spectral peakedness. 
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Figure 23. Pen-and-ink strip-chart records for selected cases with high 
values: (a) Huntington Beach record at 0100 e.s.t., 3 
September 1972; analysis at 0040 e.s.t., 3 September 1972; 
(b) Naples record at 1220 e.s.t., 9 June 1969; analysis at 
1340 e.s.t., 9 June 1969; (c) Michigan City (buoy) record 
at 2025 e.s.t., 14 September 1975; analysis at 2220 e.s.t., 
14 September 1975. 
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Table 6. Correlation between spectral-peakedness parameter and other wave parameters. 


Location Cases Correlations } Cases | Correlations 
® vs. abe Vs. Q | Tp vs. Q 


bit aun ah peo ay No. of peaks } (single-peaked cases} 


Atlantic City, N.J. ; 1,037 0.22 0.08 0.27 
Virginia Beach, Va. 0.13 0.02 QO. 
Nags Head, N.C. 0.31 0.20 0. 
Lake Worth, Fla. 0.05 0.08 0. 
Naples, Fla. 0.06 0.16 0. 
Michigan City, Ind. (buoy) 0.09 0.27 0. 
Michigan City, Ind. (pressure) 0.06 0.61 0. 
Presque Isle, Pa. (buoy) 0.10 0.14 0. 
Presque Isle, Pa. (pressure) 0.06 0.36 QO. 
Huntington Beach, Calif. 0.34 0.24 0. 
Pt. Mugu, Calif. 0.22 0.25 0. 


weakly correlated with relative water depth as evidenced by a small cor- 
relation between peak period and for single-peaked spectra (Table 6). 
This may also be a result of the definition of Qo. <A low correlation 
between significant height and Qp is also shown in the table. 


The spectral-peakedness parameter appears to be related to signif- 
icant wave height and peak spectral period in a complicated way (see 
Table 7 for Nags Head and Huntington Beach). High mean values of 
occur for significant wave heights and peak periods near the annual mean 
for most locations (Table 7). Mean values decrease fairly system- 
atically for Tp increasingly higher than the annual mean T, in most 
cases (Table 7). Q>p values have a tendency to increase with Hg for 
Tp shorter than the mean and decrease with Hs for Hg and Tp 
greater than the annual means. High values of Hg are often associ- 
ated with relatively high mean Q) for short Tp, but low mean Qp 
for long T,. The latter is mostly a result of the nonsinusoidal pro- 
file for long waves in shallow water. 


On the basis of evidence in this report and Goda (1976), it seems 
reasonable to hypothesize that high peakedness parameters are indicative 
of strong wave grouping, but low peakedness parameters in shallow water 
are not necessarily indicative of weak grouping, especially when the 
relative water depth is small. 


c. Skewness and Kurtosis of Sea-Surface Elevation Distribution 
Function. The distribution function for sea-surface elevation also pro- 
vides useful insight on shallow-water spectral statistics. Mean values 
of skewness and kurtosis (defined in eq. 6) for the selected locations 
are given in Table 8. It is evident that, except at the Great Lakes 
sites, the shallow-water distribution functions tend to have more high 
than low extremes (high skewness) and that they tend to be more sharply 
focused than the normalized Gaussian distribution (high kurtosis). These 
findings are consistent with cnoidal wave profiles shown in Figure 2. 


Cumulative distribution curves for skewness and kurtosis are shown 


in Figures 24, 25, and 26. The Atlantic coast locations have a tendency 
for higher skewness and kurtosis values than the gulf and Pacific coast 
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Table 7, Mean % values for each H, - Tp group.! 


Tp (s) Hg (cm} 
2 3 1 
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Table 8. Statistics on skewness and kurtosis of distribution of sea-surface elevations. 
Location Mean Mean Correlations 


sk i 
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Atlantic City, N.J. 
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Nags Head, N.C. 
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Naples, Fla. 

Michigan City, Ind. (buoy) 
Michigan City, Ind. (pressure) 
Presaue Isle, Pa. (buoy) 
Presque Isle, Pa. (pressure) 
Huntington Beach, Calif. 

Pt. Mugu, Calif. 
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Figure 24. 
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Cumulative probability distribution functions for skewness 
of sea-surface elevation distribution function. 
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Figure 25. Cumulative probability distribution functions for kurtosis 


of sea-surface elevation distribution function at Atlantic 
and gulf coast locations. 
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Symbol Location No. of Cases 
Michigan City(Buoy) 246 
Michigan City (Pressure) 127 
Presque Isle (Buoy) 110 
Presque Isle (Pressure) 119 
Huntington Beach 746 
Pt. Mugu 799 


Normalized Gaussian Distribution 


Kurtosis of Sea-Surface Elevation 
Distribution Function 


99.9 99 90 70 50 30 10 1 0.1 
Cumulative Percent 


Figure 26. Cumulative probability distribution functions for kurtosis 
of sea-surface elevation distribution function at Great 
Lakes and Pacific coast locations. 


locations. The Great Lakes locations and Pt. Mugu tend to have the 
lowest skewness values. The Michigan City locations and the Presque 
Isle buoy gage show a tendency for occasional very high kurtosis values 
which are associated with very low significant wave heights. Although 
several kurtosis values greater than nine occurred for the Michigan City 
locations and Pt. Mugu, the values are omitted from Figures 25 and 26. 
The Presque Isle pressure gage shows a tendency for exceptionally low 
kurtosis values. A section of pen-and-ink record for a case with high 
kurtosis (equal to 8.0) was shown in Figure 11. 


Some of the kurtosis values in Figures 25 and 26 are questionably 
high. Since kurtosis is based on the fourth power of the deviation of 
sea-surface elevations from the mean (i.e., n = 4 in eq. 6), it is more 
affected by a few elevations far from the mean than any other parameters 
considered in this report. Noise or momentary signal losses in a record 
which are comparable to the size of the largest waves will pass uncor- 
rected through the editing process described earlier. These bad data 
points usually have little effect on the spectrum but they can have a 
Significant effect on the kurtosis of the distribution of sea-surface 
elevations. 


Another factor which leads to questionably high values of kurtosis 
is the occasional tendency of the step-resistance gage to become elec- 
trically shorted at a particular elevation, often by marine growth, and 
to produce wave records in which the troughs are unrealistically flat. 
Because kurtosis is especially sensitive to these problems (known to 
have occurred intermittently in the data considered in this report), 
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the percentages associated with kurtosis values in Figures 25 and 26 
higher than about 4 should be considered somewhat inflated. 


The correlation between skewness and peak spectral period is less 
than 0.25 at all locations; for single-peaked cases it is similarly low. 
However, the correlation between skewness and significant height is 
moderatly high for the ocean and gulf coast locations (Table 8). There 
is a clear tendency for skewness to increase with significant height for 
the ocean and gulf coast locations as evidenced in a plot of mean skew- 
ness versus significant height (Fig. 27). This tendency is not evident 
for the Great Lakes locations (Fig. 28). 


It is apparent that cnoidal wave profiles (see Fig. 2) with very 
high skewness and kurtosis values are rarely observed at the selected 
gage locations. However, some of the records contain wave profiles which 
are more nearly cnoidal than sinusoidal. Figures 24, 25, and 26 indi- 
cate that cnoidal profiles appear in about 10 percent of the records or 
less. 


V. CONCLUSIONS 


The following conclusions are drawn from evidence and discussion 
presented in this report: 


1. About 60 to 70 percent of the ocean and gulf coast energy spectra 
in shallow water have one or more major peaks in addition to the dominant 
peak. 


2. At least 50 percent of the spectra in Lakes Michigan and Erie 
have one or more major secondary peaks. 


3. Major secondary spectral peaks are common for both high and low 
energy spectra. 


4. The very small sample of spectra for hurricane waves examined 
(seven spectra) does not reveal any special characteristics of hurricane- 
generated spectra. 


5. Major secondary spectral peaks are usually indicative of secondary 
wave trains. 


6. When secondary spectral peaks occur at frequencies of two or three 
times the dominant frequency, they may represent a single train of steep 
waves rather than independent wave trains, especially if the relative 
water depth for the dominant wave is small. 


7. When secondary spectral peaks occur at frequencies of 0.5 or 1.5 
times the dominant frequency and the waves are high enough for breaking 
to occur, the secondary peaks may not represent independent wave trains. 
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Figure 27. Mean skewness of sea-surface elevation distribution 
function as a function of H, for ocean and gulf 
coast locations. 
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Figure 28. Mean skewness of sea-surface elevation distribution 


functions as a function of Hg, for Great Lakes 
locations. 
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8. There is evidence that wave breaking limits the ratio H,/d at 


the shallow-water gage locations considered to less than 0.55 rather 
than the commonly used limit of 0.78. 


9. Individual shallow-water ocean spectra with comparable signifi- 
cant heights and peak periods are highly variable in shape. 


10. Mean spectra for the same significant height-peak period group 
show a surprising similarity between locations; however, there is a nota- 
ble lack of low-frequency energy in gulf and Great Lakes spectra and a 
lack of high-frequency energy in most southern Pacific coast spectra. 


11. The distribution of individual energy values about the mean 
energy for each spectral band tends to be nearly Gaussian in the vicinity 
of the dominant peak but non-Gaussian for low-energy bands away from the 
peak. 


12. For most exposed ocean locations, a spectrum is more likely to 
contain major secondary peaks when the period of the dominant peak is 
long or short than when it has intermediate values. 


13. High values of the spectral-peakedness parameter, Q» indicate 
that high waves may be occurring in groups; however, high Q, values in 
shallow water are usually associated with low-energy wave conditions. 


14. Low values of Q. generally indicate that high waves are poorly 
grouped. However, if the spectrum has secondary peaks which do not repre- 
sent independent wave trains, Q . is not a direct indicator of the extent 
of wave grouping. e 


15. Q, values tend to be higher in the Great Lakes than along the 
ocean and gulf coasts. 


16. The distribution function for sea-surface elevations in shallow 
water tends to have more high than low values far from the mean and tends 
to show less variability than the normalized Gaussian distribution. 

These features are consistent with cnoidal-type, shallow-water wave 
profiles. 


VI. SUMMARY 


Digital wave records were obtained from 11 gages at 9 shallow-water 
sites along the U.S. Atlantic, Pacific, gulf, and Great Lakes coasts. 
The records were analyzed with a fast Fourier transform procedure to 
produce energy spectra. Selected spectral parameters, including the 
number of major spectral peaks and Goda's (1970) spectral-peakedness 
parameter, and parameters of the distribution function of sea-surface 
elevations were also computed. The physical meaning of the shallow- 
water spectrum and computed parameters was discussed and illustrated 
with three cases of cnoidal wave profiles. 
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Spectra for the 24 highest energy cases from each gage were summa- 
rized individually. Many of the spectra had more than one major peak. 
Although many of the high-energy cases represent breaking wave condi- 
tions, the ratio of significant height to water depth was less than 0.55 
in all cases. This value differs significantly from the value 0.78 which 
is usually used in conjunction with depth-limited wave heights. 


All spectra were grouped according to significant height and period 
corresponding to the highest spectral peak. A mean spectrum was com- 
puted for most height-period groups. The standard deviation of spectral 
values about the mean was also included. The mean spectra show a clear 
lack of low-frequency energy at the gulf and Great Lakes sites and a 
lack of high-frequency energy at the southern California site. They also 
show evidence of systematic variations in shape as a function of signifi- 
cant height and peak period. 


Multipeaked spectra are common; 60 to 70 percent of the ocean and 
gulf coast spectra and at least 50 percent of the Great Lakes spectra 
have more than one major peak. Major spectral peaks are usually indica- 
tive of independent wave trains. 


High values of the spectral-peakedness parameter appear to be related 
to the occurrence of high waves in groups. However, low values of the 
spectral-peakedness parameter are not necessarily indicative of a lack 
of grouping of high waves in shallow water. The distribution function of 
sea-surface elevations in shallow water tends to deviate from a normalized 
Gaussian distribution in a way which is consistent with cnoidal-type wave 
profiles. 


In the course of this investigation, a computer program was developed 
to identify major peaks and valleys in an irregular signal. The program 
has many potential applications, including identification of major peaks 
in a spectrum and identification of meaningful crests and troughs in a 
time series of sea-surface elevation. 
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APPENDIX A 


INDIVIDUAL HIGH-ENERGY WAVE SPECTRAL PLOTS 
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APPENDIX B 


MEAN SPECTRA GROUPED BY HEIGHT AND PERIOD 
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Table B-1. Approximate relative depth and ratio of signif- 
icant height to depth for mean spectra grouped 
by height and period. 


Parameters based on period | Parameters based on height 
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Table B-1. Approximate relative depth and ratio of signif- 
icant height to depth for mean spectra grouped 
by height and period. --Continued 


Parameters based on period | Parameters based on height 
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Table B-1. Approximate relative depth and ratio of signif- 
icant height to depth for mean spectra grouped 
by height and period. --Continued 


Parameters = on period | Parameters based on height 
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APPENDIX C 


MEAN SPECTRA GROUPED BY HEIGHT, PERIOD, AND WATER LEVEL 


Interpretation key 


Mean Water level 
1 low 
2 mid 
3 high 
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APPENDIX D 


DESCRIPTION OF COMPUTER ROUTINE SMOOTH FOR IDENTIFYING 
MAJOR PEAKS AND VALLEYS IN AN IRREGULAR SIGNAL 


The computer routine SMOOTH is useful for deleting small, inconse- 
quential peaks and valleys from an irregular digital signal. Peaks and 
valleys that remain after application of SMOOTH represent major extrema 
which in many cases are more meaningful than small wiggles in the signal. 


SMOOTH is a very general routine which can be applied to many differ- 
ent Situations. It has been applied at CERC to the following three 
situations: 


(a) Computation of statistics on individual waves. SMOOTH was 
applied to a time series of sea-surface elevations: 


(b) Computation of statistics on tidal highs and lows. SMOOTH 
was applied to a time series of tidal elevations. 


(c) Identification of major spectral peaks. SMOOTH was applied 
to a wave energy spectrum. 


The operation of SMOOTH is most conveniently described in terms of 
its application to a time series of sea-surface elevations, although its 
other applications are analogous. The general scheme of operation con- 
sists of a check on the time difference and elevation difference between 
successive extrema. If either is less than the specified acceptable min- 
imum then one peak and one valley are deleted from the time series. 


The input to SMOOTH consists of several control parameters and an 
array (EXTIM) containing time and elevation for each extremum in the 
time series. Figure D-1 shows five extrema in a hypothetical time 
series. If the point labeled ''-1" were the first point in the time 
series, then the first 10 values in the EXTIM array would be 


t_y; N-y> to, No> ti> Ni» to, No» t3, and N35 


where t; and nz; are defined as the time and elevation associated 
with the ith point. The control parameters which must be specified 
are: 


FURST = time associated with the first point desired in the 
time series. 

ITEMS = total number of values in EXTIM array (= twice the 
number of extrema). 

CHP = minimum acceptable time difference between successive 
points (critical half period). 

HMIN = minimum acceptable elevation difference between 


successive points. 
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Cia nie) / 
“XN YY 
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Figure D-1. Hypothetical time series. 


After completion of SMOOTH, the times and elevations of major peaks 
and valleys are stored in the first ITEMS elements of EXTIM. The value 
of ITEMS has been reduced in accordance with the number of small peaks 
and valleys eliminated. Since the smoothing algorithm cannot work prop- 
erly at the end of a record, the last few points are usually accepted 
regardless of whether or not they satisfy the acceptance criteria. Thus, 
the last six elements in EXTIM may not represent major peaks or valleys 
and should either be checked or categorically eliminated. 


A complete list of subroutine SMOOTH with comments is provided in 
Figure D-2. A flow chart for the main processing loop (identified as 
loop 1520 in the listing) is provided in Figure D-3. The symbols nj 
and t; and the point numbers referred to in the flow chart are con- 
sistent with definitions in Figure D-1. Tpyp and Hyry in the flow 
chart are equivalent to CHP and HMIN in the listing. 
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SUBROUTINE SMOOTH (FURST> ITEMS) 


SUBROUTINE SMNUTH TAKES A RECORD OF PEAKS AND VALLEVS AND 
ELIMINATES YNCONSEWUENTITAL PrakKS AND VALLEYS, THE ELIMINATION 
CRITERIA ARE A MINIMUM DIFFERENCE BETWEEN PEAK AND VALLEY 
ELEVATIONS (HIN) AND A MINIMUM HORIZONTAL SPACING (BE IT 
TIMEs POSITIUONe ETC) RETWEEN PEAK AND VALLEY (CHP), 

INPUT PARAMETERS ARE DEFINED aS FOLLOWS) 

FURST & STARTING TIME 
ITEMS 3 TOTAL NUMBER OF EXTREME VaLUES CINCLUDES TIME AND 
ELEVATIUN VALUES) 
CHP 2 CRITICAL HALF PERIOD 
HVMIN © MINIMUM HEIGHT TO BE CONSIDERED 
ExTIMc(ODD) s TIME 
EXTIMCEVEN) © EXTREME ELEVATION aSSOCJATED WITH EXTIM(EVENe}) VALUE 
OF TIME 
COMMON /SMOT/ EXTIM(100) 
ig SET VALUES OF CHP AND HMIN YO BE USED 
DATA CHPe HMIN / Ve000 3,0 # 
C INITIALIZE VARIABLES AND FIND STARTING POINT FOR PROCESSING 
1503 IITEMSSITEMS 
Js} 
ITEM7TSITENSo7 
DO 1502 TEyoTITEMSo2 
TSKIJPsI 
IF CEXTIMN( I) GEeFURST) GO YO 1504 
JsJe2 
ITEMSSITEMSe2 
1502 CONTINUE 
1504 JSTARTSISKyP 


CHREKKEAEA TRATES HES SHER HHSH EA ESE REESE SKE K HE SAE KH KSEE SEK SKEKA RS SEEKS EFI ELH RESS 


MAANMARAVANANN qanaan 


¢ BEGIN MAIN PROCESSING LOOP 
DO 1520 ISYSTARTOTTEMTe2 
Ls! 
IFCISTARY.GT,JTEM7)GO TO 1520 
TFCEXTIMCTISCHP,GToEXTIMN(1¢2))G0 TH 1595 
TF CABSCEXTIMCI $1) PEXTIM(T$3)) LT HMIN)GO TO 3511 
€ IF NO TRANSFER, THIS EXTREME ACCEPTED 
EXTIM(JVSEXTIMCT) 
ExTIN¢@ Jel SExTImnc led) 


JET+2 

GO TO 1520 
Cc WHEN THE NEXT INSTRUCTION JS REACHEDe ONE HIGH AND ONE LOW WILL 
C BE DELETEN 


1591 IFCEXTIMC 191) GT ,EXTIN(1+3)IGO TO 1514 
G IF NEXT INSTRUCTION IS USEDe THIS 7$ A LOW 

IFCEXTIMC145) ,GT,EXTIM(1¢1))IGO TO 1512 

TF CEXTIMCI¢3) GT EXTIN(C1=1))GO TO 45138 


GO TO 1517 
1542 IFCEXTIMC( 143) GT EXTIM¢1+7))GO0 TO 1519 
GO 19 1516 
Cc WHEN NEXT INSTRUCTION JS USED» THIS IS 4A HIGH 


1544 IFCEXTIM( 145) STEXTI4(1+1))G0 Yo 1515 
IFCEXTIMNC 167) GT ,ExTIN(1+3))GO TO 4519 
GO TO 1516 

15145 IF CEXTIM( 73) GroEXTIM(I91))GO To 1517 
GO TO 1518 


Cc SET THE vVaLUE OF SCaSE 
1516 ICASEs{ 
GO TO 1524 
1517 ICASEs2 
Go TO 1524 
Figure D-2. List of subroutine SMOOTH. 
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1518 ICASES3 
GO TO 15214 

1519 ICASEs4 

1521 JislIeed 


Jesle7 

c DELETE ONE HIGH AND UNE LOW ACCORDING TO THE VALUE OF ICASE 
GO 10(1522515235152491525) ICASE. 

c STORAGE PLAN A 


1522 EXTI4(Je2)stxTInclea) 
ExTIM( Jeol )SExTIM(lej) 
ExTIM(J)SEXTIM(T) 
EXTIM(J+1)SExTIM(I¢)) 
EXTIM( J+e)SEXTIM( 146) 
EXTIM( J*#3)SEXTIMCI¢7) 
GU TO 1526 

c STORAGE PLAN B 

1523 EXTIM(Je2)sExTIM(Ie2) 
EXTI*(Jel)sExTIncles) 
EXTIM(J)SExXTIM( I 64) 
EXTIM(JeI)SEXTIMCIT eS) 
EXTIM(J+2)SExTIM( 146) 
EXTIM(J+3)SExTIMCI¢7) 
GG TO 1526 

c STORAGE PLAN € 

15e4 EXTIM(Je2)ysexTim( lea) 
EXTIM(JeLISENXTIMC1¢3) 
EXTIM(JISEXTIM( Je) 
EXTIM(CJ+1SEKTIMC1¢5) 
EXTI*(J+2)SExTImue leo) 
EXTIM(J+3)sEXTIMC1%7) 
GO TO 1526 

Cc STORAGE PLAN D 

$525 EXTIM(I+6)SExTIM(I42) 
ExTIM(I+7)SExTIm(le3) 
EXTIM( Jee)SsexTIn(le2) 
EXTIM( Jol) )SExTIM( lol) 
EXTIACJVSEXTIM(T) 
EXTIM( JL VSEXTIM(I91) 
ExTI4( Je2SseExTIve(Io2) 
EXTIM(C J+ SISEXTIM( 193) 
KIST+yg 
K221+40 

1526 JisJe2 
J2e5+3 
ITEMSsITEMSeg 
ISK] PsI¢8 
jstJeod 
Gc TO 1564 

1520 CONTINUE 

¢ END MAIN PROCESSING LOOP 
CHEROKEE KSEE ORES EERE ESEKERSESESES HESS SATS HETKAEESESEE SKE TH EEE CEST EE ETERS 


JeJel 
IFC Tee EU TTEMT)LES 
1530 IF(LeGTelETEMS)GO TO 1340 
JeJe! 
EXTIM(J)SEXTIMCL) 
Lele! 
GO 10 1530 
1540 CONTINUE 


ME ertnene WERE ANY DELFTIUNS IN THIS P4&SS REPROCESS aLL REMAINING 


HIGnS AND LO#S TO Make FURTHER DELETIONS IF NEEDED. 
IFC(ITITEMS.GT,1TEMSIGO TO 1503 
RETURN 


ND 5 
Figure D-2. List of subroutine SMOOTH .--Continued 
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Point Exit 
0 is within 3 
points of 
record en 


processing 


loop 


4 P number points so 
Point 0 is iS P 


that point 
accepted P 


2 becomes point 0 


ICASE = 1 ICASE = 4 


Remove points Remove points Remove points 
1 and 2 0 and 1 -1 and 0 


Renumber points so 


Remove points. 
2 and 3 


that point 


4 becomes point 0 


Figure D-3. Flow chart for main processing loop in SMOOTH (loop 1520). 
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